Optimal squeezing and entanglement from noisy Gaussian operations 
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We investigate the creation of squeezing via operations subject to noise and losses and ask for 
the optimal use of such devices when supplemented by noiseless passive operations. Both single and 
repeated uses of the device are optimized analytically and it is proven that in the latter case the 
squeezing converges exponentially fast to its asymptotic optimum, which we determine explicitly. 
For the case of multiple iterations we show that the optimum can be achieved with fixed intermediate 
passive operations. Finally, we relate the results to the generation of entanglement and derive the 
maximal two-mode entanglement achievable within the considered scenario. 



Squeezed states are a valuable resource for different 
fields of physics. They can increase the resolution of pre- 
cision measurements, as exploited in gravitational wave 
detectors Q, improve spectroscopic sensitivity Q, and 
enhance signal-to-noise ratios 0, e.g., in optical commu- 
nication. Moreover, squeezing acts as a basic building 
block for the generation of continuous variable entan- 
glement 0, which in turn is a cornerstone for quan- 
tum information purposes. Unfortunately, squeezing 
is an expensive resource as well: squeezed states are 
hard to create and the involved operations are subject 
to losses and noise inevitably restricting the attainable 
amount of squeezing. On the other hand passive opera- 
tions, in quantum optical setups implemented by beam- 
splitters and phase shifters, can often be performed at low 
cost and they are — compared to the squeezers — relatively 
noiseless. 

This work is devoted to the question, how can we ex- 
ploit a given noisy squeezing device in an optimal way 
when supplemented by arbitrary noiseless passive opera- 
tions. We derive the optimal strategy for single and re- 
peated use of the squeezing device, calculate the achiev- 
able squeezing and relate it to the maximal attainable 
amount of entanglement. To this end we will use a black 
box model for the physical squeezing device. This will 
give us the possibility to derive optimality results which 
are equally applicable to a wide range of physical imple- 
mentations. 

The argumentation will make use of the covariance 
matrix formalism, which was mainly developed in the 
context of continuous variable states having a Gaussian 
Wigner distribution — so called Gaussian states The 
latter naturally appear in quantum optical settings (the 
field of a light mode) as well as in atomic ensembles (col- 
lective spin degrees) and ion traps (vibrational modes). 
We restrict ourselves to the natural class of Gaussian op- 
erations, i.e., operations preserving the Gaussian char- 
acter of a state 0, 0] • This includes all time evolutions 
governed by operators quadratic in bosonic creation and 
annihilation operators. All the presented results hold for 
an arbitrary number of modes and although it might be 
reasonable to think in terms of Gaussian states, we do 
not have to restrict the input states to be Gaussian. 



Preliminaries. — We will begin with introducing the 
notation and recalling some basic results [E El Ll El- 
Consider a system of N bosonic modes with respective 
canonical operators (Qi, Pi, . . . ,Qn, P/v) = R- These 
are related to the annihilation operators via aj = (Qj + 
iPj)/y/2 and satisfy the canonical commutation rela- 
tions [Rk, Ri] = *Cfe;l governed by the symplectic matrix 
°i qJ ■ The displacement din phase space 

and the covariance matrix (CM) 7 corresponding to a 
state p are then given by 



(1; 



tr[pRi] and 7^ = tr [p{(Pj - dj), (Rj 



dj)}+] 



where {•,•}+ denotes the anti-commutator. 

While for coherent states 7 = 1, a state is called 
squeezed if its uncertainty in some direction in phase 
space is below the uncertainty of the coherent state, i.e., 
if s(7) = A m in(7) < 1, where 5(7) is the squeezing of 
7 measured by its smallest eigenvalue (9J. Note that by 
this definition, more squeezing means a smaller s. As 
the squeezing is independent of the displacement rf, we 
omit it for the remaining part of the paper. 

Let us now focus on Gaussian operations. Unitary 
Gaussian operations are precisely those realizable by 
quadratic Hamiltonians, so that they naturally appear 
in many physical systems. In phase space they act as 
symplectic operations S £ Sp(2iV) on the covariance 
matrix 7 1— > S T jS. Symplectic operations preserve the 
canonical commutation relations and are thus given by 
the group Sp(2JV) = {Se m. 2Nx2N . S T aS = a y An 
important subgroup is given by the group of orthogonal 
symplectic transformations K(2N) = 0(2N) H Sp(2N) 
[9|. Physically, these correspond to passive operations, 
which can, in quantum optical set ups, be implemented 
by beam-splitters and phase shifters [IOJ . Obviously, pas- 
sive transformations cannot change the squeezing, since 
elements from K(2iV) preserve the spectrum and in par- 
ticular the smallest eigenvalue of the CM. 

We will now introduce the model we use to describe the 
squeezing device. In general, a noisy operation £ can be 
regarded as a noiseless map on a larger system including 
the environment, which is discarded afterwards. If the 
overall time-evolution is governed by a quadratic Hamil- 
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FIG. 1: Various scenarios for the optimization of squeezing: 
a) single iteration case: from a given input ji n> we want to 
generate as much squeezing as possible by properly choos- 
ing K and using the noisy device £ ; b) multiple iteration 
case: the device can be applied repeatedly, and we have to 
determine the Ki's such as to maximize the finally obtained 
squeezing; c) circular setup: with identical Ki = K. 



tonian and the environment is in a Gaussian state (e.g., a 
thermal reservoir), it can be shown that these operations 
are exactly the ones which act as 

7 i-> £ (7) = X T jX + Y , X T iaX + Y >ia (1) 

with X,Y e R 2Nx2N 6]. Here, the equation on the 
right hand side ensures complete positivity, i.e., guaran- 
tees that the operation is physically reasonable. While 
the X Tr yX part of £ represents a joint rotation and dis- 
tortion of the input 7, the Y contribution is a noise term 
which may consist of quantum as well as classical noise. 
In the following, we will consider squeezing devices of the 
type in Eq. as these are the ones which naturally ap- 
pear in many experiments. We will, however, show at the 
end of the paper that the results even hold for arbitrary 
Gaussian operations (which may include measurements 
and conditional operations). 

The question we are going to investigate is the fol- 
lowing: given some noisy device £ and the set of passive 
operations K(2TV), how can we generate squeeezing as ef- 
ficiently as possible from a given input state? Naturally, 
this general question can be asked in several specific ways. 
First of all, one might ask how much squeezing can be 
generated by a single application of the device given a 
certain initial state, as in Fig. 

The much more interesting question, of course, relates 
to an iterative scenario, i.e., how can we generate squeez- 
ing as efficiently as possible by repeated application of £ 
with passive operations Ki inbetween (Fig. ^3)? In this 
scenario we may either allow for different Ki or choose 
them identically as it is for instance the case in a ring 
cavity setup (Fig. [TJ;) . 



Single iteration. — In order to prepare for the more 
complicated scenarios, let us first have a look at the case 
of a single iteration (Fig.QJi); starting from a given input 
CM 7. This is the basic building block for all the iter- 
ative protocols. We now use a formal trick in order to 
facilitate the derivation: we split the CM into two parts 

(Fig. UK), 



7 = si + N 



(2) 



where s = s(j) is the squeezing of 7. The first part can 
be regarded as a "coherent kernel" of 7. It may have sub- 
Heisenberg variance (if s < 1) and it is invariant under 
passive operations. The second part is a "noise term" 
N which ensures that 7 is a physical state. As s is the 
smallest eigenvalue of 7, N > has a null space. 

Let us now see what happens if we rotate 7 = si + TV 
by some passive operation K and then send it through 
£: the coherent kernel is invariant under K , and thus 



7 



£(K T jK) = £ (si) + X T K T NKX , 



i.e., the action of the squeezing device on the "coher- 
ent kernel" plus the "noise" part which has been ro- 
tated and squeezed by X. Note that the first part no 
more depends on the choice of the passive operation; fur- 
thermore, as TV > and thus X T K T NKX > 0, the 
smallest eigenvalue of £ (si) gives a bound to the squeez- 
ing of the output. In the following, we show that this 
bound can be achieved. Therefore, let so = s(£ (si)) = 
A m ; n [sX T X + Y] be the squeezing obtained from the in- 
put si with corresponding eigenvector \<jj) [llj. On the 
one hand, s' = 3(7') > sq, and on the other hand, 



s' < (</>]£ {K T 1 K)\<j>) = <0|£(sl) 
= a<H 



f X T K T NKX\ 
X T K T NKX 1 6) . 



Recall that by definition @, TV has a null eigenvector 
which we denote by \v), N\v) — 0. By choosing K such 
that K(X\(f>)) oc \v) (which can be always done with 
passive K, cf. |3]), the second term vanishes, and we 
indeed find that s' = so = s(£ (si)). 

The proof is also illustrated in Fig. [5] Choosing K ap- 
propriately ensures that the noise TV does not contribute 
in the most squeezed direction. Note that for a given £ 
it is now straight forward to derive an optimal K and 
by exploiting the results of 01 to decompose it into an 
array of beam-splitters and phase shifters. 

For a single iteration of £ this shows that the opti- 
mally achievable squeezing s' at the output is given by 
the squeezing obtained from the non-physical input si: 



s' = f £ {s) :=X mhl [sX T X + Y] 



(3) 



It is highly interesting to note that, therefore, the op- 
timal value of the final squeezing does only depend on 
the initial squeezing (and on the properties of £ ), but on 
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FIG. 2: Illustration of the single iteration optimality proof. 
Left: the input state is decomposed into a coherent kernel si 
and into some extra noise N with a null eigenvector \v) . Right: 
after the application of the operation £, the coherent kernel is 
squeezed in the direction \<f>). By choice of K, one can achieve 
K(X\4>)) oc \u), i.e., the noise N leads to no contribution in 
the most squeezed direction. 

no other property of the initial state, irrespective of the 
number of modes considered. It can be easily checked 
that the respective function fg in Eq. yj is concave, 
monotonously increasing, and fg(Q) > (Fig. |3J. 

Multiple iterations. — The fact that the optimal output 
squeezing does only depend on the input squeezing im- 
mediately implies that in the case of multiple iterations 
the passive operations Ki can be optimized successively 
in order to obtain the global optimum. This is a remark- 
able result, as in general problems of this kind require 
optimization over all parameters (i.e., over all Ki) simul- 
taneously. Graphically, the squeezing in each iteration 
step moves along a zig-zag line between fg(s) and the 
identity, as shown in Fig.|3 The optimal output squeez- 
ing Sqo (for number of iterations — + oo) is determined by 
fs {soo) = $oo, Sqo > 0. By inserting the definition JSJ of 
fg and solving for we obtain 

-1 

The convergence of s(7i) to the optimal value is expo- 
nentially fast and bounded from above and below by the 
slope of fe at Sqo and the slope from (soo, fe(s oo)) to the 
starting point (so, fe(so)), respectively. Note that for 
s(lm) < Soo, however, s(7 ou t) > s(7in), so that squeezing 
is destroyed by applying the operation £ . 

In realistic physical scenarios, it might be difficult to 
tune the passive operations independently and it is more 
likely that the same physical device will be passed again 
and again, e.g., in a ring cavity (cf. Fig. ^i), and thus 
only one fixed passive operation K = Ki Vi can be im- 
plemented. In the following, we demonstrate that this 
is already sufficient to reach the optimal squeezing s^. 
The proper K is the one which preserves the squeezing 
at the optimality point, corresponding to a zig-zag line 
along the tangent of fg at Soo. The convergence is thus 
still exponentially fast. 

In order to see how this works, consider the non- 
physical output 

7 = £( Soc l) = Soa X T X + Y 
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FIG. 3: The function s' = fs(s) describes the optimal output 
squeezing s' which can be obtained from input squeezing s 
by an operation £ . The straight line is the identity, and the 
gray zig-zag line describes the evolution of / when using the 
optimal strategy for multiple iterations, cf. Fig. 03. 

obtained from the input Sool. By the properties of Soo, 
it is clear that A m i n (7) = with a corresponding nor- 
malized eigenvector \ipoc)- Now choose Koo such that 

|Voo) <X Xco^oo) ■ (4) 

This is exactly the K which preserves the optimality 
point, as IV'oo) is the null eigenvector of 7 — SqqI. For any 
initial state 7 with 5(7) > Soo, choose the decomposition 
70 = Sool + P , where P > and A min (-P) + s^ = s = 
s(7o). After one iteration 71 = £ (K^joKoo), we have 

+ {i)oo\X T K^PK^X^oo) 

^ J 
V 

a 

As we will show in a moment, a < 1, and for the squeez- 
ing s n = s("f n ) after n iterations it holds by recursion that 
s n < s 00 + a n (ip 00 \P\ip 00 ), which converges exponentially 
to Soo- From 

Soo(^oo\X T X\lp oc ) + ( , oo |y|V'oo) = Soo 

it follows that 

a = (ipoolX A|-0oo) = 

Soo 

which is positive and strictly smaller than one as long as 
(4><x\Y\ipoo) > 0, which is generically the case |l2j . 

Example. — Let us now consider an example which illus- 
trates how the representation of the operation in Eq. Q 
is related to the master equation p — i[p, H] + C[p] of a 
system. From it, one obtains a master equation for the 
evolution of the covariance matrix, 

7 = A7 + jA T + N . (5) 

For the case of photon losses to a vacuum reservoir, 
C[p] = v(2apa^ — a) ap— pa^ a) , one obtains A = 1o~H — v\ 
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FIG. 4: Optimization results for a physical device St which 
is given by applying the master equation for a time t, with 
H = §a f a - ±aV + h.c. and v = 0.1 (cf. text). The dotted 
line shows the squeezing obtained by simple applying £ t to a 
coherent state, while the solid line gives the asymptotically 
attainable optimum as derived in the paper. The dashed line 
shows the rotation angle which leads to the optimal value. 



and N — 2v\, where H is the Hamiltonian matrix, i.e., 
H = (Q,P)H(Q,P)1. By integration, one finds that ap- 
plying JSJ for a time t leads to a map £ t : 7 1— » X T jX + Y 



with X = e~ ut exp[-2iJcrt] and 



Y = X 1 



2v 



e 2uT e -2*HT e 2H*T dT 



X 



At this point, the results derived in the paper can be 
applied. Fig. ^ shows a one mode example with TL = 
ja^a — ^a'a> + h.c. and the noise level v = 0.1. Note 
that additional passive operations enhance the attainable 
squeezing although the noise N is rotationally invariant. 

Entanglement generation. — The optimality result on 
squeezing generation has direct implications for optimal 
entanglement generation. Indeed, it has been shown that 
the amount of entanglement which can be generated by 
passive operations starting from two squeezed Gaussian 
input states only depends on their squeezing irrespective 
of the number of modes [13j: for two modes, e.g., this 
is done by sending the states onto a beam splitter af- 
ter rotating them into orthogonal directions. Using this 
result, we can immediately determine how much entan- 
glement we can generate from Gaussian inputs using a 
black box £ supplemented by passive operations, namely 
^V.opt = l°g( s opt); where s op t is the maximal squeezing 
generated by £ and the entanglement is measured by the 
logarithmic negativity |l4J |. Moreover, by combining the 
results it is straightforward to explicitly derive the opti- 
mal entangling protocol for any given Gaussian device. 
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Appendix. — In the following, we show that the results 
obtained for channels of the type also hold for the 
most general type of Gaussian channels which may in- 
clude measurements and postprocessing. Channels of 
this type appear, e.g., in the creation of spin squeezing 
using quantum nondcmolition measurements with feed- 
back [l5( • The most general memoryless operation on N 
modes can described by a 2N mode covariance matrix T 
via the Jamiolkowski isomorphism Q as 



£ (7) = A - C(B + 7 r 1 C T , where T = 



A C 
C T B 



and r is the partial transpose of T. Again, s = s( 7 ), \ip) is 
the eigenvector corresponding to the smallest eigenvalue 
of £ (si), and we need to show that K can be choosen 
such that s(£ (si)) = s{£{K 1 K T )), i.e., that (ip\C(B + 
sl)" 1 ^^) = (i)\C(B + K 1 K T )- 1 C T \ip). This means 
that for |x) = C T \xp), P = (B + si) > 0, and N = 
7 — si > with a null eigenspace we have to find a K 
such that 



(xlP^-OP + z^vi^r 1 !*) 



(x\p 



-1/2 



1- 1 + P- 1/2 KNK T P- 1/2 



(6) 

p- 1/2 \x) 



(*) 



vanishes. Since (*) has the same null eigenspace as 
P- 1 / 2 KNK T P- 1 I 2 , the expression © can be indeed 
made zero by choosing K such that \v) cx K T P^ 1 \x) 
(where N\v) = 0). This proves that for any Gaussian op- 
eration the optimal output squeezing can be computed 
on si, thus generalizing the results of the paper. 
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